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Quark Confinement in the Deconfined Phase 
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In cylindrical volumes with C-periodic boundary conditions in the long direction, static quarks are confined even 
in the gluon plasma phase due to the presence of interfaces separating the three distinct high-temperature phases. 
An effective "string tension" is computed analytically using a dilute gas of interfaces. At T c , the deconfined- 
deconfined interfaces are completely wet by the confined phase and the high-temperature "string tension" turns 
into the usual string tension below T c . Finite size formulae are derived, which allow to extract interface and string 
tensions from the expectation value of a single Polyakov loop. A cluster algorithm is built for the 3-d three-state 
Potts model and an improved estimator for the Polyakov loop is constructed, based on the number of clusters 
wrapping around the C-periodic direction of the cylinder. 
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1. INTRODUCTION 

It has long been established that the Z(3) 
global symmetry of the pure SU (3) gauge theory 
gets spontaneously broken at high temperatures 
jlj, leading to three distinct bulk phases, sepa- 
rated by interfaces. In this talk, we will present 
an unusual quark confining mechanism due only 
to the presence of these interfaces in a cylinder, 
relating the free energy of a single static quark to 
these domain walls. We believe that this is also 
evidence that the interfaces observed in numeri- 
cal simulations correspond to physical gluonic do- 
main walls in Minkowski space-time. 

2. INTERFACE GAS 

Pure Euclidean SU(3) gauge theory at finite 
temperature T = 1//3 is described by the ac- 
tion S[A„] = f£dtfd 3 x (l/2e 2 )Tr F^F^ and 
is periodic in the Euclidean time direction. The 
Polyakov loop <&(:?) is constructed from the Eu- 
clidean time component of the gauge field A±(x, t) 
and under gauge transformations g(x, t) which 
are also periodic in Euclidean time, both the ac- 
tion and the Polyakov loop are unchanged. How- 
ever, if the gauge transformations differ by a cen- 
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ter element i.e. g(x, t + (3) = g{x,t)z,z e Z('S), 
the action is invariant but the Polyakov loop 
changes into Q'(x) — &(x)z. The expectation 
value ($} = exp(— (3F) measures the free energy, 
F, of a static quark. In the confined phase, F di- 
verges and ($) vanishes, while in the deconfined 
phase F is finite and ($} ^ 0. Hence, the Z(2>) 
center symmetry is spontaneously broken at high 
temperatures. 

With spontaneous symmetry breaking, the spa- 
tial boundary conditions and the manner of the 
infinite volume limit are important. We consider 
a spatial volume of size L x x Ly x L z . If we apply 
periodic boundary conditions, ($) vanishes even 
in the deconfined phase, due to Gauss' law in a 
periodic volume [|| — topologically, we cannot 
have a single quark in a periodic box because its 
center electric flux cannot go to infinity, it must 
end in an anti-quark. 

We apply C-periodic boundary conditions in 
the z-direction only [||. When a C-periodic field 
is translated by L z , it is replaced by its charge 
conjugate. For example, for C-periodic gluons, 
Afj,(x + L z e z , t) = A^(x,t)*, where * means com- 
plex conjugate. Physically, we can have a sin- 
gle quark in a volume, partnered with an anti- 
quark on the other side of the C-periodic bound- 
ary. Imposing C-periodicity explicitly breaks the 
Z(3) symmetry [Q, but this symmetry breaking 
disappears in the infinite volume limit. With C- 
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periodic boundary conditions, ($} is always non- 
zero in a finite volume. 

The original motivation for applying C- 
periodicity in the long direction of a cylinder was 
to extract the string tension a from numerical 
simulations measuring ($) in the confined phase. 
A single static quark sits in the cylinder, con- 
nected by a tube of gluons to its anti-quark part- 
ner on the other side of the C-periodic boundary. 
Then we have ($) = S exp(— (3aL z ) and so the 
free energy of a single quark in the cylinder is 

F = log So +<tL z . 

There are several advantages to using this method 
in a numerical study. Firstly, it is much easier to 
measure the Polyakov loop than a Wilson loop or 
some other correlator. Secondly, in a periodic box 
of size L, the Wilson loop cannot be larger than 
L/2 in its spatial extent, whereas this technique 
allows us to exploit the entire volume. 

We can also investigate (<3>) in the deconfined 
phase at temperatures T > T c , where the three 
distinct deconfined phases coexist ||. They are 
distinguished by different values of (<E>) — one 
expectation value is = ($ ,0), which is ro- 
tated by Z(3) transformations to give $ ( - 2 - ) and 
$( 3 ). A typical configuration in a cylinder con- 
sists of several bulk phases, aligned along the z- 
direction, separated by deconfined-deconfined in- 
terfaces. The interfaces cost free energy F pro- 
portional to their area A = L x L y , such that the 
interface tension is given by add — F/A. The ex- 
pectation value of the Polyakov loop can be cal- 
culated from a dilute gas of interfaces ||. The 
interface expansion of the partition function can 
be viewed as 
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The first term has no interfaces and thus the 
whole cylinder is filled with deconfined phase d± 
only. An entire volume filled with either phase d 2 
or e?3 would not satisfy the boundary conditions. 
The second and third terms have one interface 
separating phases d 2 and d$. Here, C-periodic 
boundary conditions exclude phase d\. For each 



configuration, we integrate the Boltzmann weight 
over all possible locations of the interface, then we 
sum the Boltzmann weights for all configurations 
with all possible numbers of interfaces. Summing 
the interface expansion to all orders, we obtain 
Z = exp(-(3f d AL z + 2jexp(-f3a d dA)L z ). Here, 
fd is the bulk deconfined free energy density and 
7 is a factor resulting from capillary wave fluc- 
tuations of the interfaces. Note that in three di- 
mensions, 7 is to leading order independent of 
the area A 0. To calculate ($), we simply in- 
clude the Polyakov loop expectation value of each 
configuration with the Boltzmann weight in the 
interface expansion. Summing to all orders, we 
obtain ($) = $ exp(— 37 exp(— )3addA)L z ), from 
which we calculate that the free energy of a single 
static quark in a C-periodic cylinder is given by 

1 37 
F = -~log$ + -£exp(-/3a dd A)L z . 
P P 

This result is counter intuitive. Although we are 
in the deconfined phase, the quark's free energy 
diverges in the limit L z — > 00, as long as the cross 
section A of the cylinder remains fixed. This is 
the behavior one typically associates with confine- 
ment. In fact, a' = (3^/(3) exp(— flotddA) plays 
the role of the "string tension" , even though there 
is no physical string connecting the quark to its 
anti-quark partner on the other side of C-periodic 
boundary. Because the deconfined-deconfined in- 
terfaces can lead to the divergence of a quark's 
free energy and have physically observable conse- 
quences, we believe that they are more than just 
Euclidean field configurations. 

Because the phase transition is of first order 
||, as we approach T c , the confined phase can 
coexist with the three deconfined phases, so we 
can also have confined-deconfined interfaces with 
an interface tension a c d- At T c , there are two 
possibilities [p|Jl0|l. If we have add = 2a c d, a 
deconfined-deconfined interface always splits into 
two confined-deconfined interfaces, separated by 
a film of confined phase — this is called com- 
plete wetting. If add < 2oi c d, both deconfined- 
deconfined and confined-deconfined interfaces are 
stable — this is called incomplete wetting. Nu- 
merical simulations indicate that the gluon sys- 
tem undergoes complete wetting ]ll|,|q|. In that 
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the interface expansion of the partition func- 



c 


+ 


di 


c 


di 


c 


+ 


di 


c 


d* 



} + 



The sum over i extends over the three deconfincd 
phases and d* denotes the charge-conjugate of di. 
For example, d* = d\ and d\ = d$. Note that 
due to complete wetting, one always has an even 
number of interfaces. As before, we sum the ex- 
pansion to all orders to calculate the partition 
function and (4>), obtaining 



$o cxp(f3xL z ) + E exp(-/3(a: + a)L z ) 
2 cosh[/3L z ^Jx 2 + (3(5 2 //3 2 )exp(-2/3a C(i A)] ' 



Here, S is the factor characterizing capillary wave 
fluctuations of the confined-deconfined interfaces 
and x — |(/ c — fd)A measures the bulk free en- 
ergy difference between confined and deconfined 
phases. From this expression, we can extract an 
effective "string tension" a' which, as we lower 
the temperature further into the confined regime, 
reduces to the standard string tension a in the 
large A limit, as expected. 

We can also consider the incomplete wet- 
ting scenario by including configurations with 
deconfined-deconfined interfaces at T w T c . 
The effective "string tension" we extract again 
matches a as we move into the confined regime. 

3. IMPROVED ESTIMATOR 

Using the analytically derived finite size for- 
mulae, we can determine the interface tensions 
add and a c d from numerical simulations. Be- 
cause simulating SU (3) gauge theory is computa- 
tionally intensive, we select a simpler model, the 
3-d three-state Potts model. Its action is given 
byS= where 3> x e Z(3). 

This model has a Z(3) symmetry and a first or- 
der phase transition between three distinct or- 
dered (deconfined) phases and one disordered 
(confined) phase. Here, ($) corresponds to the 
average spin. We have developed a single cluster 
algorithm with an improved estimator for mea- 
suring ($} in a cylinder with C-periodicity im- 



posed in the long direction. As we grow the clus- 
ter, we keep account of how many times the clus- 
ter crosses the C-periodic boundary. If the cluster 
wraps around the cylinder in the long direction 
an odd number of times, we count a 1. Oth- 
erwise, we have a non-wrapping cluster and we 
count a 0. This is because a non-wrapping clus- 
ter is not aware of the C-periodicity and so can 
flip to any of the three deconfined phases, making 
on average a zero contribution to ($), whereas a 
wrapping cluster can only be in deconfined phase 
d\ satisfying C-periodicity. Because we have an 
improved estimator, we can accurately determine 
very small quantities. Work on simulations to ex- 
tract the interface tensions is in progress. 
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